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DETECTOR 



.AMPLIFIED PHOTODtODE 
(SLOW-RESPONDING) 



Samples such as thin polymeric fiJms (20) are analyzed using optically induced phonons by excitation of the sample 
using radiation (32, 34) preferably absorbed by the sample and probe radiation (48, 76) preferably not absorbed by the sample, 
that is diffracted from the surface of the sample. The pulse width of ihe probe is preferably on the order of the detectable 
diffraction signal so that the phonon decay from each excitation pulse can be detected (58, 72) and analyzed (64, 62, 37). The 
technique is applicable to various samples by inducing a ripple morphology on the sample surface and detection of light 
diffracted substantially from surface ripple. 
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MEASUREMENT OF MATERIAL PROPERTIES 
WITH OPTICALLY INDUCED PHONONS 

Field of the Invention 
5 This invention relates to the meastirement of 

material properties using optically induced phonons. 

Background of the Invention 
Laser Induced Phonons (LIPS) are produced by time- 
coincident laser pulses intersecting inside a sample, 

10 setting up an optical interference pattern, i.e., 

alternating intensity peaKs and nulls. Energy deposited 
into the system via optical absorption or stimulated 
Brillouin scattering results in the launching of 
counterpropagating ultrasonic waves (phonons) whose 

15 wavelength and orientation match the interference pattern 
geometry. The mechanism by which LXPS ultrasonic waves 
are generated depends upon whether the sample is 
optically absorbing or transparent at the excitation 
wavelength. If the excitation pulses cure absorbed e.g. 

20 into high-lying vibronic levels, rapid radiationless 
relaxation and local heating at the interference maxima 
(the transient grating peaks) occurs. Thermal expansion 
then dirives material in phase away from the grating peaks 
and toward the grating nulls, setting up 

25 counterpropagating waves. 

In samples which are transparent at the excitation 
wavelength, optical energy is coupled directly into the 
sample's acoustic field via stimulated Brillouin 
scattering. This process takes advantage of the inherent 

30 spectral line width in 100-picosecond (psec) excitation 
pulses. Higher-frequency photons from each pulse are 
annihilated to create lower- frequency photons in the 
opposite pulse and phonons of the difference frequency 
and wave vector in the medium, Counterpropagating waves 

35 (a standing wave) are thus produced. 
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In either case the acoustic wave propagation, 
which continues after the excitation pulses leave the 
sample, causes time-dependent, spatially periodic 
variations in the material density, and since the 
5 saunple's optical properties (real and Imaginary parts of 
the index of refraction) are density-dependent, the 
irradiated region of the sample acts as a Bragg 
diffraction grating. This propagation of the optically 
excited ultrasonic waves can be optically monitored by 

10 time-dependent Bragg diffraction of a variably delayed 
probe laser pulse. 

Summary of the Invention 
An object of the invention is to make LIPS 
(alternatively referred to herein as Impulsive Stimulated 

15 Thermal Scattering or ISTS) measurements by reflectance 
of a probe beam from the sample surface without regard to 
the absorbance of the probe beam by the sample. Another 
object is to make LIPS measurements without regard to 
sample type or thiclcness. A particular object is the 

20 analysis of thin, polymer samples by LIPS-ref lectance 
with a non-absorbing probe laser. 

The phenomenon which allows realization of these 
objects is believed to be the inducement of a physical 
surface morphology or 'ripple' in response to laser 

25 excitation. It has been discovered that the surface 

ripple phenomenon can be selectively induced based on the 
nature of the sample, e.g. sample stiffness and thickness 
and the angle of incidence and wavelengths of the 
excitation radiation and selectively detected or analyzed 

30 by control or monitoring of parameters such as angle of 
incidence, polarization and intensity of a probe beeun 
reflected from the ripple morphology. The discovery 
enables advantages in terms of the samples that can be 
analyzed, including ultrathin polymer films, in terms of 

35 the rate at which data can be obtained, such as obtaining 



BNSDOCID: <WO 9301476A1_L> 



SUBSTITUTE SHEET 



wo 93/01476 



PCr/US92/05679 



- 3 - 

the full transient signal induced by each excitation 
laser shot, and in terms of the type and accuracy of data 
that can be obtained, including more accurate subtraction 
of probe pulse contribution from the signal and the 
5 collection of new data such as acoustic dampings rates 
(attenuation), thermal conductivity and film 
delaminat ion . 

In a particular example, excitation of a sample 
using optically absorbing wavelengths can be used to 

10 induce the physical morphology, preferably in the form of 
a periodic grating on the surface of a sample. The 
grating morphology can be interrogated by reflectance of 
a probe laser from the surface without regard to the 
probe laser wavelength and in particular, allowing the 

15 use of a probe laser not substantially absorbed by the 
sample. The sample may be of various types. A 
particular example is the analysis of thin polymer films. 
The films, which may be UV^absorbing and visible- 
transmitting, may be excited by UV pulses and 

20 interrogated by reflectance with visible pulses. For 

very thin samples, e.g., 500, 30, 10 or even l/xm or less, 
this operational mode is particularly advantageous since, 
with a non-absorbing probe laser, sample heating which 
can easily damage thin samples, is reduced. It is also 

25 an advantage that samples can be interchanged without the 
need for selecting new probe laser wavelengths. 
Alternatively, in the case of thin samples, e.g., polymer 
films, acoustic waves may also be excited through 
impulsive Brillouin scattering in which case the 

30 excitation pulses need not be absorbed by the sample. 

Another object is to meOce measurements quickly, 
with fewer laser shots, to make it possible to observe 
changing sample properties with high time resolution. 
Employing techniques of the invention, LIPS measurements 

35 may be made in real time. An excitation pulse having a 
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short pulse width compared with the phonon oscillation 
period may be employed in combination with a high power 
probe laser having a relatively long pulse width such 
that a substantial time-portion, in many cases all, of a 
5 sample's transient time dependent response from each 
excitation pulse can be measured. While a high power 
probe may be particularly useful in some circumstances, 
alternatively, it has been discovered that signal can be 
obtained using a relatively low power probe, e.g. about 1 

10 watt, provided by a CW source such as an argon ion laser. 
In either case, the probe signal from successive 
excitation pulses may be signal averaged to increase 
signal to noise. The detection electronics are selected 
so they have sufficiently fast response times to permit 

15 time resolution of the phonon oscillations* Operation in 
these modes is, again, particularly advantageous for 
thin, polymer samples since the number of excitation 
pulses can be reduced and thus the number of heating 
cycles, which can potentially damage thin, fragile 

20 samples can also be reduced. For a probe ii^ich has a 
short width compared to the lifetime of the transient 
response, a temporal delay between excitation and probe 
pulses may be varied to interrogate the sample at 
different time segment stages in its response evolution. 

25 The term "polymers" as used herein refers to 

molecules composed of sequences of repeating monomer 
units connected by covalent bonds. A particular polymer 
is generally not made up of a single molecular species. 
Rather, polymers are mixtures of macromolecules with 

30 similar structures and molecular weights that exhibit 

some average characteristic properties. The polymers may 
be homopolymers or copolymers, linear or branched, cross- 
linked or bonded by electrostatic interaction such as 
hydrogen bonds. The polymers may be elastomers, 

35 thermoplastics and may include plasticizers, stabilizers. 
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anti-oxidants and lubricants. Polymers may be crystalline 
in the sense that long segments of linear polymer chains 
are oriented in a regular manner with respect to one 
another. Such crystalline regions of a polymer are 
5 referred to as crystallites. Amorphous, noncrystalline 
regions generally lie between the crystallites and 
constitute defects in the crystalline structure. 

A "sample" such as of a polymer may include 
additives , e.g., chromophores incorporated within the 
10 polymer mass (typically not covalently bonded) for 
enhancing absorbance of the sample at a particular 
wavelength such as a desired excitation or probe 
wavelength. A "pure sample" as used herein contains no 
such additive. 

15 The term "absorption", as used herein, includes 

absorption of electromagnetic radiation by ground or 
excited states of the sample molecule. Particular 
aspects of the invention use LIPS in the reflectance mode 
where the absorbance of the probe is less than the 

20 absorbance of the excitation. Particularly the probe 
absorbance is 50%, 10% or even 1% or less than the 
excitation absorbance. In some cases, no absorbance of 
the probe beeun by the sample can be detected. Nolar 
absorptivities (£) at the probe wavelength may be less 

25 than 1000 or 100. Absorption may be less than about 10% 
of the maximum absorption (which may be relative to 
complete absorbance) in a wavelength region (e.g., UV, 
visible, infrared), typically less than 1%. The 
excitation wavelength is selected such that the 

30 absorbance is sufficient to give rise to a periodic 
morphology on the sample surface, detectable as a 
diffraction signal. 

In a particular aspect, the invention features an 
apparatus for measuring the properties of a sample of 

35 material. The apparatus includes a first, excitation 
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source for producing excitation radiation adapted to 
impinge upon the sample of material. The excitation 
radiation is comprised of pulsed radiation composed of at 
least two component pulses which interfere within the 
5 sample* The excitation radiation is sufficient to induce 
a transient phonon in the material which gives rise to a 
transient, time dependent periodic ripple morphology on a 
surface of the sample. A second, probe source is 
provided for producing probe radiation arranged to 

10 reflect from the periodic ripple morphology on the 

surface of the sample to form a diffraction signal. A 
detector detects the diffraction signal from the probe 
source radiation reflected from the surface, and an 
analyzer selectively analyzes the diffraction signal 

15 formed by the transient ripple morphology. 

In another particular aspect, a method for 
measi2ring the properties of a sample. The method 
includes impinging a pulse of excitation radiation upon 
the sample. The excitation radiation is composed of at 

20 least two component pulses which interfere within the 
sample and are selected to induce a transient phonon in 
the sample which gives rise to a transient, time 
dependent ripple morphology on a surface of the sample. 
Probe radiation is reflected from the periodic ripple 

25 morphology on the surface of the sample to form a 

diffraction signal. The diffraction signal from the 
probe source radiation reflected from the ripple 
morphology is selectively detected. 

In another particular aspect, the invention 

30 features an apparatus for measuring the properties of a 
thin sample of polymeric material. The apparatus 
includes a first, excitation laser source for producing a 
pulse of radiation adapted to impinge upon the sample. 
The excitation source is a pulsed source composed of at 

35 least two component pulses which interfere within the 
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thin film. The xcitation radiation is adapted to induce 
a transient phonon in the material. A second, probe 
laser source is provided for producing radiation. The 
probe radiation is of selected wavelength not 
5 substantially absorbed by the sample and arranged to 
reflect from the surface of the sample to form a 
diffraction signal. A detector detects the diffraction 
signal. 

The features of the above aspects can be combined. 

10 Apparatus features can be used in method inventions. In 
addition, in various aspects the invention may include 
one or more of the following features. The radiation 
from the probe source is absorbed about 50%, 10% or 1% or 
less than the radiation from the excitation source or 

15 there is no detectable or substantial absorbance by the 
sample at the probe wavelength. The excitation source of 
radiation is in the ultraviolet. The probe source is in 
the visible. The sample is a thin sample, e.g., on the 
order of about 500, 30 or 10 fim or less, and the 

20 excitation radiation is absorbed by the sample. The 
sample is a thin polymeric film. The film is a free- 
standing film or disposed on a support. The probe 
radiation is diffracted from the surface of the sample. 
The probe beam and detector are for detecting a 

25 substantial time-portion of the time^dependent 

diffraction signal, induced by each excitation pulse such 
as the entire detectable time-dependent diffraction 
induced by each excitation pulse. The probe radiation is 
of a selected pulse width and the detector is adapted to 

30 detect the time dependent diffraction for the duration of 
the probe pulse. The excitation pulse width is on the 
order of psec duration and the probe pulse width is on 
the order of nanosecond (nsec) duration (e.g. generally 
greater than 10 nsec, e.g. around 100 ns) . The detector 

35 has a time resolution on the order of 1 nsec. A signal 
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averaging means is provided for signal averaging the 
diffracted radiation from multiple excitation pulses. The 
probe pulse has a peak power output (power during the 
laser on-time) of about 1000 watts or greater, such as 
5 around 10,000 watts. The probe laser is a Q-switched YAG 
laser. The excitation pulse is generated from a pulsed 
laser and the probe pulse can be generated from a CW 
laser having for example a laser power output of around l 
watt, such as a gated argon ion laser. The polymer 

10 sample is a pure polymer sample. 

Various aspects may also include one or more of 
the following features. The diffraction signal is 
analyzed to selectively analyze diffraction from surface 
ripple. The analyzer includes a polarizer for 

15 determining change of polarization of the probe beeun 

after diffraction to selectively analyze diffraction from 
the surface morphology. The analyzer may also be for 
analyzing the signal as a function of wavevector to 
selectively analyze diffraction signal from the ripple 

20 morphology. The angle of incidence of the probe beam 
and/ or the excitation beam or the wavelength of the 
excitation beam may be varied to optimize diffraction 
from the ripple morphology. The diffraction signal is 
produced by reflection of the probe beam from the back 

25 surface of a sample, opposite the radiation sources or 
the front surface or both. The system is adapted for 
determining adhesion and/or thermal diffusion of the 
polymer sample on a substrate surface from the 
diffraction signal. 

30 The invention has many advantages and 

applications. It is demonstrated below that ISTS 
experiments with real-time data acquisition rates are 
possible and can be applied to studying the pseudo- 
Ray leigh acoustic modes that propagate in thin film 

35 coatings. Real-time data acquisition can be crucial in 
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avoiding optical damage to polymer films. The diffracted 
ISTS signal eirises predominantly from corrugation or 
ripple displacements at the film surface and the 
film/ substrate ripple due to the propagating pseudo- 
5 Rayleigh modes. The relative contribution of each 
acoustic mode to the ISTS signal depends on the 
efficiency with which each mode is excited by the 
excitation pulses and on the efficiency with which each 
mode diffracts the probe beam due to surface and 

10 interface ripple. A formalism is developed to quantify 
both of these factors. Using this formalism a method for 
extracting the elastic constants of the film by analyzing 
the dispersion of the pseudo-Rayleigh modes is described 
and applied to the polyimide/silicon system. 

15 The techniques have wide applications in 

nondestructive, real-time characterization of thin films. 
For many applications, such as monitoring thin-film 
fabrication and cure or determination of the spatial 
uniformity of coatings, it is not necessary to determine 

20 elastic moduli but merely to monitor changes in acoustic 
frequency. Finally, characterization of thermal 
diffusion rates will present a similar range of 
applications in polymeric, diamond, and other thin films. 
Other types of materials can be studied, for example, 

25 single crystal Fe films and amorphous, plasma deposited 
carbon films. 

The invention includes methods of use of the 
techniques described. Still other objects, features and 
advantages follow. 

30 Detailed Description 

We first briefly describe the drawings. 
Drawings 

Figs. 1 to Ic illustrate schematically, the 
surface ripple morphology phenomenon included in the 
35 present invention; 
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Fig. 2 is a schematic of a measurement apparatus 
according to the invention; Fig. 2a is an absorhance 
spectrum of a polymeric sample; Fig. 2b is a plot of the 
intensity verses time, raw data, from the experimental 
5 setup in Fig. 2; and Fig. 2c is a plot of the processed 
data; 

Fig. 3 is a schematic of an alternative 
measurement apparatus according to the invention; Fig. 3a 
is a plot of the intensity verses time, raw data, from 

10 the experimental setup in Fig. 3; Fig. 3b is a plot of 
the probe pulse intensity; Fig. 3c is a plot of the 
processed data from Fig. 3a and Fig. 3d is a plot of 
processed data taken under the conditions of Fig. 3a, but 
from a single excitation shot; 

15 Fig. 4 is a plot of intensity versus time raw 

data, and Fig. 4a is a plot of the processed data from an 
experimental set up similar to Fig. 3, with an acoustic 
wavevector of 0.60 faoT^; 

Fig. 5 is a plot of intensity versus time, raw 

20 data, from an experimental set up using a CW probe laser 
while Fig. 5a is a plot similar to Fig. 5 but from a 
single excitation shot; 

Fig. 6 is a plot of intensity versus time, raw 
data, from an experimental set up using a CW probe laser 

25 to obtain thermal data; 

Fig. 7 illustrates the system geometry used in the 
theoretical calculations, wherein the film fills the 
space between y=0 and y=-h and the substrate fills the 
area from y=0 to y=+H and the acoustic wavevector set up 

30 by the pump beams is in the z direction; 

Fig. 8 is a plot of pseudo-Ray leigh mode 
dispersion generated using substrate elastic parameters 
v^a = 8945 m/s, v^g = 5341 m/s, Pg = 2.33 g/cm*^ and film 
elastic parameters v^f = 1300 m/s, v^g = 700 m/s, pg = 

35 1.45 g/cm-^, with the acoustic wavevector given by q, the 
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film thickness is h, and v^g and Vj^^ are the substrate and 
film Rayleigh surface velocities' respectively; 

Fig. 9 is a schematic illustration of the beam 
paths for the four components of first order ripple 
5 diffraction that contribute to ISTS signal, wherein parts 
A and B describe the different paths for the two orders 
of diffraction observed experimentally; 

Fig. 10 illustrates material displacements, 
generated using the optimized (for the system studied) 

10 elastic constants, for the first eight pseudo-Rayleigh 
modes on a l^m film coating an infinite substrate at a gh 
value of 2.5, with the scale factors (A) for each mode 
determined such that the largest component of 
displacement in the mode which is driven most efficiently 

15 is equal to one and the ripple amplitude (R) for each 

mode given relative to the maximum displacement amplitude 
for that mode; 

Fig. 11 is a plot of the film surface ripple 
amplitude versus qh and velocity for a given ISTS heat 

20 input using the optimized (for the system studied) 
elastic constants wherein each symbol represents a 
different mode dispersion curve to facilitate comparison 
with Fig. 14 and the units of the z axis are arbitrary; 

Fig. 12 is a plot of the film-substrate interface 

25 ripple amplitude versus qh and velocity using the same 
parameters as Fig. 11 and the units of the z axis are the 
same as those given in Fig. 11; 

Fig. 13 is a plot of lattice distortion of the 
first two modes at gh=0.8, which, upon comparison with 

30 Fig. 11 indicates the existence of a crossover; 

Fig. 14 is experimental data (symbols) and 
theoretical pseudo-Rayleigh mode dispersion curves (solid 
lines) for pyralin films on silicon substrates with 
various film thicknesses (h) in microns; the elastic 

35 parameters used in generating the dispersion curves were 
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op-binized to best fit the data that fall on the lowest 
two velocity curves; 

Fig. 15 is a plot of dispersion curves for the 
first 30 normal modes for a 3^m film on a 330/im substrate * 
5 using the optimized elastic constants for the 
pyralin/silicon system; the existence of modes is 
illustrated with velocities above v^g which was the cut- 
off value for the pseudo-Ray leigh modes calculated for an 
infinite substrate as in Figs. (6) and (12) ; 

10 Fig. 16 is a dispersion curve plot of velocity 

versus wavevector for a tightly bound film- substrate 
system, while Fig. 16a is a similar plot for a 
freestanding film; 

Fig. 17 is a series of dispersion ctirve plots of 

15 velocity versus wavevector; 

Fig. 18 illustrates the system geometry for an 
analysis of thermal diffusion. 
General Description 

Referring to Figs. 1-lc, a sample 2, typically a 

20 thin T^, e.g. 5/tm thick or less, polymer thin film is 
provided having a front exposed surface 3 and a back 
surface 5 which is disposed on a substrate 6. (It will 
be understood that the film could be, as well, a free 
standing film, without substrate 6.) As shown in Fig. 

25 la, the sample is excited by impinging a pair of 

excitation beams 7,9 from sources E^^, Ej. The wavelength 
of radiation from and £2 ^® selected such that it is 
absorbed by the sample. The beams 7,9 interfere within 
an overlap region 11 within the bulk of the sample, 

30 causing rapid local heating in a periodic fashion. Upon 
such excitation, a periodic ripple- like morphology 13 is 
generated at the surface 3 . The morphology has a peak to 
trough height, H, that is at a maximum soon after 
absorption of the excitation wavelengths and then 

35 oscillates and decays with time. In particular 
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experiments, as discussed below, maximum heights H, are 
estimated on the order of 5-10 nanometer (nm) . As shown 
in Fig. lb, the morphology 13 may be interrogated by a 
probe laser pulse 15, from a source P^, by diffraction of 
5 the pulse in the reflection mode from the surface 3. The 
diffracted beam 15' is detected by a detector D, 
positioned on the same side of the sample as the source 
Pj^. The angle of incidence of the probe, e^^, and angle of 
diffraction e^, may be varied and selected based on the 

10 spacing of the periodic morphology. The source P^^ is 

preferably selected such that the wavelength of the probe 
pulse is not substantially absorbed by the sample. As 
illustrated, a component 15" of the beam 15 may be 
transmitted through the sample. The pulse 15 is 

15 typically selected such that its duration is long 

compared to the excitation pulse and the lifetime of the 
periodic morphology. The pulse may be from a pulsed 
probe laser of selected pulse width and energy or a 
continuous wave laser, for example. The detector is 

20 adapted to receive the signal pulse 15', which is the 
diffracted beam component from the sample surface, for 
analyzing the characteristics of the morphology over the 
course of its existence from each excitation. As 
illustrated in Fig. Ic, after analysis, the sample 2 

25 returns to its original state. 

The. excitation wavelength is preferably in the 
ultraviolet (about 100-4 oonm) , but may be, e.g. , in the 
visible (about 410-800nm) , infrared (about 0.1-50/x) or 
beyond. The probe is preferably in the visible, but 

30 maybe in other regions, e.g. infrared or ultraviolet, 

depending on sample absorbance. The excitation and probe 
radiation may be in the same wavelength region. An 
advantage of the invention is that the sample may be 
probed with a laser wavelength independent of the sample 

35 absorption; greater flexibility is provided in the choice 
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Of the probe laser. In particular, it is advantageous to 
employ a probe laser emitting visible light which can be 
easily aligned in manufacturing applications. As further 
discussed in the examples, analysis times can be reduced 
5 by employing a probe pulse width as long as the transient 
decay of the induced phonons, so that the entire 
morphology decay is detected from each excitation. 

A particular advantage is the analysis of thin 
films I especially thin polymer films, e.g. lO/zm thickness 

10 or less, or even l/im thiclcness or less. Typical polymer 
samples include visible wavelength transmitting polymers 
such as polyethylene, polyurethane and polyimides. 
Polymer films which may be analyzed include molecular 
films such as Langmuir-Blodgett films or biological 

15 films. The films analyzed may be supported on a 

substrate or may be freestanding films. In the case of 
films on a substrate, the support provides both 
mechanical rigidity and can reduce sample deunage when the 
support acts as a heat sink. Sample supports include 

20 silicon or glass. Applications for the system include 
emalysis of polymer protective coatings on silicon 
substrates and various applications to polymers in 
polymer processing, including in-situ processing during 
e.g. extrusion. In the latter case, the analysis may be 

25 carried out at various points in the extruder to 

determine the progress of the operation and to modify the 
extrusion in response to the data obtained. Other 
polymer applications include polymer curing and polymer 
loading (analysis of properties under mechanical load] • 

30 The properties that can be measured by analysis of 

the raw data include acoustic speeds, determined from the 
frequency of surface oscillations; attenuation rates 
determined from decay of the signal; thermal expansion 
and thermal diffusion rates; thermal expansion 

35 coefficients and heat capacities can be determined from 
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signal intensities cases where stsuidard samples are run; 
compressibility and elastic modulus may be determined 
from acoustic speed and sample thickness frequency 
dependencies. Calculation of many properties from the 
5 raw data is known in the art, and is further discussed in 
Journal of Chemical Physics 94, 7677 (1991) and Farnell 
et al. , "Elastic Wave Propagation in Thin Layers", 
Physical Acoustics, (W-P. Mason, R.N. Thurston Eds.) 
935FF, Academic Press, NY, 1972. A particular method for 
10 measuring thermal diffusion is discussed below. 

Calculation of excitation efficiencies is also discussed 
below. 

Example 1 

Short Probe Pulse 

15 Referring to Fig. 2, an apparatus 2 for measuring 

the properties of a thin film of polymeric material 20, 
such as a thin (about 2.2 micron) film of Pyralxn^ 2555 
(a polyimide polymer available from E.I. DuPont de 
Nemour, Wilmington, DE; other Pyralin materials may be 

20 also available for study, e.g. PI 2525, PI 2545, PI 2611) 
on a silicon substrate 22, includes a first laser 24 such 
as a YAG laser that produces a beam at a wavelength of 
about 1.064 micron and an energy of about 1 millijoule. 
The pulse width is on the order of picoseconds, in this 

25 example 100 psec. The beam is modified by crystals, 

26,28 for frequency doubling and tripling to produce an 
excitation beam wavelength in the ultraviolet at about 
355 nanometers. Another preferred wavelength is about 
266 nanometers. The beam from the laser is split by a 

30 beam splitter 30 into component beams 32, 34 which are 
" directed by reflecting optics 36, 38 toward a focusing 

element 40 (about 30 cm focal length cylindrical focus) 
^ which focuses component beams 32, 34 such that they 

overlap in time and location within the sample 20. 

35 Typically, the spot size on the sample is selected to 
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avoid excessive hea'ting and damage tiherefrom, e.g., "the 
spot size nay be 100 microns to l millimeter with an 
energy of about 10 /tJ. in this example, the spot size 
was about 100 fim high and 2 mm wide (cylindrical focus) • 
5 The probe laser 42 is a YAG laser for probing the surface 
morphology by reflectance. The output beam from the 
laser 42 passes through a Pockel's cell 44 to isolate a 
particular pulse with an energy of, for example, about 50 
microjoules. The probe beam is passed through a doubling 

10 crystal 46 to produce a probe pulse 48 having a 

wavelength generally not absorbed by the sample 20 such 
as a probe pulse in the visible range, e.g. 532 nm, as in 
this example. The probe pulse 48 is directed to the 
sample by optical reflectors 50, 52 and focused by 

15 focusing element 54 (about 20 cm focal length) such that 
the beam impinges upon the sample surface at angle 
about 45^ from normal to the unperturbed sample surface. 
Diffracted beams 56, eurising from the morphology at the 
sample surface, may be detected at angle Bj^ about 45^, by 

20 a detector 58 such as an euaplified photodiode. The 
pulses from the excitation laser 24 and probe laser 42 
are electronically delayed to probe the surface 
morphology as a function of time. A delay control 60 
times the sequential pulses of the lasers. Typically, 

25 the excitation laser pulses at about 1 kilohertz. The 
probe laser pulses at delay times (from the excitation 
pulse) varying from the nanosecond range to 40 
microseconds or more, depending on the properties of the 
sample to be determined. The probe pulse width is short 

30 relative to the period of the phonon oscillations so that 
each probe pulse thus represents a single measurement of 
the instantaneous diffraction signal. The signal 
detected by the detector 58 is filtered via a lock-in 
amplifier 62 and passed to an analysis means 64, e.g., a 

35 microprocessor, which also receives the pulse delay 
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information and which converts the raw data via Fourier 
transform to a spectrum of intensity versus frequency. 

Referring to Fig. 2a an absorbance spectrum of a 
thin film of Pyralin<» is shown. At the excitation 
5 wavelength of 355 nm, the material is highly absorbing, 
s At the probe wavelength of 532 nm, on the other hand, 

radiation is not substantially absorbed. 

Referring to Fig. 2b, the raw LIPS data for the 
thin Pyralin^ sample is shown as a plot of the intensity 
10 versus time, being the delay between the excitation beam 
pulse and the probe beam pulse. In Fig. 2c, the 
processed data is shown after compilation with the Fast 
Fourier Transform algorithm (FFT) . The analysis time for 
this measurement was about 10 minutes using about 150,000 
15 excitation and probe pulses. The data may h& analyzed to 
yield material properties as described above. 

Example 2 

Long Probe Pulse 

Referring now to Fig. 3, in an alternative mode, a 

20 probe laser 70 is adapted to provide a probe beam having 
a relatively long pulse width such that a substantial 
time-portion of the transient grating can be measured 
from each excitation pulse. The probe pulse is typically 
of longer duration than the excitation probe pulse and 

25 may be as long or longer than the detectable diffraction 
signal from the sample surface. In this example, like 
numbers refer to the same elements as descriloed in 
Example 1, Fig. 2. The probe laser 70 may be a YAG laser 
operated in the Q-switch mode with energy of about 3 

30 mini joules in a pulse width of about 300 nanoseconds. 
(The peak power, therefore, was on the order of 10,000 
watts. Lower power operation, ^^.g., down to around 1000 
" watts is believed to be highly practical. Power 

requirements may be dependent on signal strength which 

3 5 may vary with the sample.) In this case, the probe laser 
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beam 76 impinges upon the sample surface simultaneously 
with the excitation beams and the full time dependent 
response from each excitation pulse is recorded. (If 
signal persists after 300 ns, the probe pulse can be 
5 delayed as discussed in Example 1.) The detector 72 is a 
fast response (370 psec) amplified photodiode (model ARX- 
SA, Antel Optronics, Inc. , 3325B Mainway, Burlington, 
Ontario, CANADA) • The signal is received by a transient 
digitizer 74 vith a l GHz bandwidth (digital signal 

10 analyzer model 602, Tektronix, Inc., P.O. Box 500, 
Beaverton, OR 97077) . The i GHz bandwidth and fast 
photodetector provide about 1 nsec time resolution, so 
that acoustic waves of up to 1 GHz frequency can be time 
resolved. The detector and electronics are generally 

15 selected to be as fast as the sample response of 
interest. A polarizer (thin film) is preferably 
positioned between lens 54 and the sample. 

Referring to Fig. 3a a plot of raw data the 
transient signal from a single excitation pulse is shown. 

20 An advantage of operation in this mode is that, rather 
than probe a single time point of the transient 
morphology as in the embodiment of Fig. 2, the 
diffraction signal from the full duration of the 
transient morphology produced by each excitation pulse is 

25 detected. The time required for signal detection is 

therefore considerably reduced. Fig. 3a also illustrates 
the probe pulse profile (which can be subtracted from the 
raw data prior to processing) . In Fig. 3b, the processed 
data (without subtraction of the probe pulse profile) is 

30 shown. This data was collected in about 15 seconds, 

using about 500 excitation shots, the signals from which 
were collected and signal averaged. Pyralin films of 
various thicknesses have been studied, e.g. from about 
0.92 to 5.81 micron (spin coated on a 330 micron thick 
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silicon wafer with thickness determined within ±0.05 
micron by mechanical stylus prof ilometer) • 

Referring particularly to Fig* 3a, the 
oscillations in the data are due to acoustic oscillations 
5 initiated by the crossed excitation pulses. The signal 
intensity depends on the induced time-dependent material 
displacements (which are described in more detail in the 
sections that follow) , weighted by the probe pulse 
intensity envelop which is shown in Fig. 3b. in Fig. 3a, 

10 the first part (-85 ns) of the probe pulse arrive at the 
sample before the excitation pulse, so there is no signal 
other than that due to parasitically scattered probe 
light. The excitation pulses heat the sample and 
initiate acoustic oscillations which are observed in the 

15 data. There is also a nonoscillatory component of the 
signal due to steady— state thermal expansion in the 
heated regions of the sample. This component finally 
decays as thermal diffusion washes away the spatially 
periodic variation in sample temperature. Thermal decay 

20 is not observed in these experiments because at these 

wavevectors, the timescale for thermal diffusion is much 
longer than the length of the Q-switched probe pulse. 
Instead, the signal in Fig. 3a disappears as the probe 
pulse ends. (As discussed below, with a larger temporal 

25 probe pulse, such as a cw probe pulse, longer sample 
excitation signals may be detected.) 

From the data in Fig. 3a, the frequency of the 
surface acoustic wave can be extracted. In the Fourier 
transform of the data shown in Fig. 3c, the acoustic 

30 response is dominated by one mode with a frequency of 113 
+ 2 MHz. The dc component of the transform (not shown) 
is 30 times larger than the acoustic peak. In Fig. 3d 
the Fourier transform of diffracted signal is shown taken 
under identical conditions as the data in Fig. 3a, but 

35 with only one laser shot, i.e., no averaging. While this 
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spectrum is noisier than the spectrum of the averaged 
data, the position of the main peak is identical. (The 
extra noise seen in Fig. 3d is believed to be due to 
self-nodelocking within the one Q-switched probe pulse.) 
5 The data clearly illustrates the ability to perform true 
one-laser shot experiments. Best results can be obtained 
with a smooth probe pulse temporal profile. 

In Figs. 4 and 4a, ISTS raw data and its Foxirier 
transform, are shown respectively for the same sample at 

10 a larger scattering wavevector. It can clearly be seen 
that there are at least two surface acoustic modes which 
contribute to the signal. 

To insure that the signal was due to the polyimide 
film, control experiments were performed on uncoated 

15 silicon wafers. No signal from the silicon surface could 
be detected. 

Further, it is possible to distinguish between 
scattering due to modulation of the dielectric constant 
through elasto-optic coupling and diffraction due to 

20 surface and interface ripple through analysis of the 
polarization properties of the diffracted light and the 
acoustic wavevector dependence of the diffracted signal 
intensity. Pseudo-Ray leigh waves have both transverse 
and longitudinal character, and therefore give rise to 

25 modulation of both diagonal and off -diagonal elements of 
the dielectric tensor of a thin film through elasto-optic 
coupling. If elasto-optic coupling is a significant 
source of diffraction, then a component of the diffracted 
probe beam should be polarized at £m angle 90*> relative 

3 0 to the incident probe beam (depolarized or VH 

scattering) . The magnitude of this component relative to 
the non-depolarized diffracted component depends on the 
square of the ratio of off -diagonal to the diagonal 
elasto-optic coupling coefficients. By contrast, for a 
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coplanar scattering geometry surface ripple does not 
rotate the polarization of the diffracted light. 

In addition, if elasto-optic coupling were the 
dominant diffraction mechanism, the signal intensity 
5 would be independent of scattering wavevector q since the 
acoustic strain induced through bulk ISTS excitation is 
g- independent. By contrast, surface ripple arises from 
components of the acoustic displacement rather than 
strain and should thus vary as g"^. Signal from surface 

10 ripple depends on the square of the displacement, and as 
detailed below the signal intensity due to surface ripple 
decreases approximately as g~^. 

In experiments performed above, no depolarized 
component of the diffracted signal was detected. It was 

15 also observed that the diffracted signal intensity 

decreases sharply as the q is increased. This indicates 
that the diffracted signals are due predominantly to 
surface and interface ripple induced by the pseudo- 
Ray leigh waves. 

20 Example 3 

CW Probe 

A further embodiment employs a continuous wave 
probe pulse from, e.g., a CW argon laser generating a CW 
beam of known duration in place of the YA6 laser in 

25 Example 2. The probe laser wavelength may be in the 
visible, 514 nm, and probe pulse width controlled by an 
£0 modulator (in place of the doubling crystal shown in 
Fig. 3). (No electronic delay is necessary for the CW ' 
experiments.) This embodiment, similar to that of 

30 Example 2, allows the transient response from each 
excitation shot to be completely detected. In this 
method the probe beam has a time- independent profile, 
thus enabling reduction in noise contribution when 
subtracting the probe beam profile from the detected 

3 5 signal. 

SUBSTITUTE SHEET 



BNSDOCID: <:WO 9301476A1 I > 



wo 93/01476 



PCT/L'S92/05679 



- 22 - 

Referring to Figs. 5 and 5a, raw data from an 
unsupported polyimide film using a YAG laser excitation 
so\irce at X=26e nm, 100 psec (a quadrupling crystal is 
used in place of tripling crystal 28) and a gated (500 
5 ^sec) small frame argon ion CW probe (about 50 100 

micron beam size, with single line power of about 1 watt) 
is illustrated for 200 signal averaged excitation shots 
and a single excitation shot, respectively. Despite the 
relatively low intensity, about 50 to 100 times lower 

10 than the Q-switched probe used in the experiments above, 
good data was obtained. Care should be taken to 
efficiently gate the beam to avoid overexposure which can 
lead to damage of the sample- In these experiments, a Q- 
polarizer was placed before the £0 gate. 

15 An advantage of the gated CW probe is that thermal 

process which occur on timescales larger than acoustic 
processes can be efficiently monitored. Likewise, the 
long temporal widths available with the CW probe allows 
the entirety of the sample acoustic profile to be 

20 detected. In addition, by including frequency dependent 
loss modulus, the experiments can account for acoustic 
dampings. Referring to Fig. 6, a thermal decay curve of 
the polyimide sample over a microsecond time scale is 
shown. Similar data can be obtained from a single 

25 excitation shot. 

Theoretical Treatment and Parameter Optimization 

A discussion of the features of the data obtained 
and a theoretical analysis that leads to optimization of 
the experiment for detecting signal predominately from 
30 the induced surface ripple, follows. 

In ISTS, two ultrashort laser excitation pulses 
are crossed temporally and spatially at the sample. An 
optical interference pattern characterized by scattering 
wavevector q is formed, where q is the difference between 
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the vavevectors of -the crossed pulses. The scattering 
wavevector magnitude g associated with each excitation 
angle 6^ can be calculated according to: 




5 where A, is the wavelength of the excitation light. 
Experimentally, referring particularly to the data 
obtained in the experiments performed as discussed with 
respect to Fig. 3 et seq., data can be recorded as a 
function of the excitation angle 6 for all the Pyralin/Si 

10 seunple thicknesses. The angles used in this study are 
0.48*, 0.68», 0.92«, 1.07*», 1.17«», 1.62* 1.67», 

1.93*, 2.37«», 3.40", 3.92«», and 5-85*. (These 
angles may be meastired mechanically using a calibrated 
rotation stage and are accurate to +-0.05*. ) The 

15 corresponding scattering wavevectors according to (a) are 
0.15, 0.21, 0.22, 0.28, 0.33, 0.36, 0.50, 0.52, 0.55, 

0. 60, 0.73, 1.05, 1.21, and 1.81 inverse microns. These 
are accurate to +0.01 inverse microns. The wavevectors 
were chose to lie along a crystal axis of the silicon 

20 substrate. However, experiments performed with a 

wavevector magnitude of 0.52Atm"^ on a 5.81/Lm Pyralin film 
sample showed no dependence on wavevector direction. 

1. Surface Ripple Amplitude 

An order of magnitude estimate can be made for the 
25 surface ripple generated in an ISTS experiment by 

calculating the magnitude of the temperature grating set 
up by the excitation laser pulses and then relating this 
to a displacement amplitude using the linear thermal 
expansion coefficient. The temperature change aT for 
30 light impinging on the surface (in the y--h plane) in the 
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y direction and setting up a grating of wavevector q in 
the z direction can be written as in eqs (2-3). 



AT(y.2) = Ac-^>[l+c'^^J 



(2) 



1 



A=^Pa(l-R)rECe-0> (3) 



5 Here, 2.303C is the material absorption coefficient, Ig 
is the total energy per unit area of the excitation laser 
pulses, p is the mass density, C is the heat capacity per 
unit mass, R is the reflectivity for the material/air 
interface and is the fraction of absorbed light that 

10 is converted to heat* An upper limit for the average 
value ot the temperature increase over the thicJcness of 
the film can be estimated with (2) by setting =1 and 
R=0. For these experiments, Is-l 6 mJ/cm^ and C»1.3^m'^. 
Using values of p = Ig/cm^ and C 2.3J/(g*K) which are 

15 typical for polymers leads to an average temperature 
increase of approximately 8 K for a 4/Jim thick film. 

This temperature change can be related to the 
change in the y-direction length of a volume element 
S'Ly{y,z) using the linear thermal expansion coefficient 

20 as in eq. (4) . 

5Ly(y.z) = aTAT(y.2)dy (4) 

For the general case, the total surface ripple Hfi.p(z) can 
be found by integrating this expression over the sample 
thickness. For a thin film of thickness h attached to a 
25 nonabsorbing substrate as in these experiments the 

integration extends only over the thickness of the film. 
Performing this integration leads to eq. (5) • 
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Rrip W = ^ Pa (1 - R)Ie(i - c'^'Wl + C^^] ( 5 ) 



Using a typical polymeric value of 0^= 80X10"^ in/(in.K)^® 
for the linear thermal expansion coefficient yields an 
upper-limit estimate of 0.005 microns for the surface 
5 corrugation amplitude excited by ISTS. This estimation 
assumes that neighboring volume elements slide by one 
another, and does not consider the effects of stress that 
will tend to reduce the corrugation. For this reason, 
expression (5) does not show the expected 1/g dependence. 

10 II. Theorv 

In order to extract the elastic constants of the 
thin film from the pseudo-Rayleigh mode frequencies 
measured at different scattering wavevector magnitudes q, 
the dispersion of uCq) for the various modes is to be 

15 understood. In fact, the frequency depends only on the 
product qh, where h is the film thickness, so that the 
results of measurements on films of different thickness 
can be compared. As mentioned earlier, the ISTS 
excitation efficiencies and the probe diffraction 

20 efficiencies for each of the modes sure to be determined. 
The ISTS signal expected from a supported film is derived 
below treating the excitation then the probing process. 
An isotropic theory is used since no effects of the 
elastic anisotropy of silicon or the film were observed. 

25 A. ISTS Excitation of the System 

Starting with the equations of thermoelasticity 
for an isotropic, elastic, and homogeneous medium: 



8^u 

C44V^u + (cu-C44)V(V.u) = yVT-hp-^ 
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(7) 



10 



15 



In these expressions u is a vector describing the 
material displacements, T is the fluctuation in 
temperatiire relative to the equilibrium temperature, p is 
the equilibrium density, and C44 and c^^ are the elastic 
constants which are related to the bulk longitudinal and 
transverse acoustic velocities according to Vj^=(C2i/p) 
and Vt-(C44/p) , y is a constant related to the elastic 
constants and the coefficient of linear volume expansion 

by Y = ( 30^^-4044 ja^, k is the thezmal conductivity, 
^=yTo/K and Cg is the constant strain specific heat per 
unit mass. Finally, Q represents the absorbed heat per 
unit time per unit volume derived from the excitation 
laser pulses. 

If u is expressed in terms of longitudinal and 
transverse (?) potentials. 



then equations (6) and (7) can be rewritten as follows. 



Vy = 0 



(8) 




(9) 




ciii? pen 



T 



(10) 




C44 



(11) 
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It is assumed that longitudinal compressions do not cause 
substantial changes in temperature (r^-'O in eg. (7)). 

Note that for a bulk system, the last two 
equations describe uncoupled longitudinal and transverse 
5 acoustic modes respectively. Since the laser heating 
enters through eq. (9) which is only directly coupled to 
eq. (10) , ISTS only excites longitudinal acoustic modes 
in bulk isotropic systems. However, this simple picture 
breaks down for a thin film, since boundary conditions 

10 must be satisfied. In particular, for a thin film in 
intimate contact with a substrate, the material 
displacements and the normal components of the stress 
must be continuous across the interface and the normal 
components of the stress at any free surface must vanish. 

15 These boundary conditions lead to coupling between the 
longitudinal and transverse potentials so that every 
surface acoustic mode has both longitudinal and 
transverse character. 

The film/ substrate system geometry is depicted in 

20 Fig. 7. The substrate fills the region from y=0 to y=+H 
and the film fills the region from y-O to y=-h. It is 
assumed that the ISTS excitation beams form an infinite 
uniform grating pattern with wave vector <i in the z 
direction which is independent of x and damped along y 

25 due to optical absorption. In this analysis, H is 

considered infinite so that the substrate fills the whole 
y>0 half space. (The corrections that occur when the 
finite size of the substrate is explicitly accounted for 
is discussed below.) 

30 The first step in solving equations (9) -(11) is to 

determine the temperature distribution set up by the 
excitation pulses. The dynamics of thermal diffusion 
take place on a much longer timescale than the acoustic 
dynamics. Thus, for analyzing the excitation of the 

35 film-substrate system, k can be set to zero in eq. (9) • 

SUBSTITUTE SHEET 



BNSDOCID: <WO_9301476A1J_> 



wo 93/01476 



PCT/US92/05679 



- 28 



In ISTSr the excitation laser pulse duration is short 
compared to the acoustic oscillation period. Here it is 
assumed that conversion of optical energy to heat through 
molecular electronic and vibrational relaxation also 
5 occurs on a fast time scale relative to the acoustic 
oscillation period so that the time dependence of Q in 
eg. (9) can be approximated by a delta-function at zero 
time. More gradual thermal izat ion will reduce all the 
acoustic amplitudes but should not affect the relative 
10 amplitudes of the different modes substantially. 

Equation (9) is solved for the temperature distribution 
set up by the laser excitation beams to yield eg. (12). 



Here, e(t) is the Heaviside step function which turns on 
15 at t=0, and A and C were introduced in egs (2-3) • since 
the film absorbs strongly and the silicon substrate 
reflects the 355nm excitation light used in these 
experiments, T(y,z,t) is only nonzero in the film (for - 
h<y<0) . 

20 With this functional form for the temperature 

distribution, equations (lO) and (11) can be solved using 
transform techniques. Since the temperature distribution 
is independent of x, all derivatives with respect to x 
are neglected. From the geometry of the system and the 

25 nature of the temperature distribution given in (12), a 
Laplace transform for the time variable and exponential 
Fourier transform for the z coordinate are the natural 
choices. (Notation is such that S is the Laplace 
variable conjugate to t and k is the Fourier variable 

30 conjugate to z.) The transformed solutions for the 
potentials <p and ? are then found to be. 
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A^(k.s)cxp|y^k2 + 
5(k) 



♦(>c.y.$) = 



ps l-s^ + 



C^^ -k'vtf +C''vtf-s''J 



+ B^(k.s)exp -y k2 + 



-h < y < 0 







f I 2 ^ 




+ D^(k,$)exp 









H>y>0 



(13) 



.2 



Ay(k.s)xexp| yjk^ ^^^''•^^^"P 



-h < y < 0 



f Hz $ 

Cy(k.s)xcxply^k'' + - 



2 



+ Dy(k.s)xcxpf-yjk2+^ 



H>y>0 



(14) 



where k is tihe Fourier 'transform variable conjugate ^o z 
and s is the Laplace variable conjugate to time. In 
5 these expressions^ ^Lf , ^Tf , ^Ls, and ^Ts are the bulk 
longitudinal and transverse velocities in the film and 
substrate respectively. In taking the Laplace transform, 
the initial conditions that there is no displacement or 
motion at t=0 when the laser pulses first arrive at the 
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saunple have been used. Note that choosing the only 
nonzero component of T to lie along the x axis allows for 
trivial satisfaction of the gauge v»i^=0 used here. 
Finally, all but the third term in the expression for the 
5 longitudinal potential in the film are homogeneous 

solutions to the equations of motion. This third term is 
the particular solution of equation (10) and can be 
thought of as driving the motion of the film and giving 
rise to a response consisting of a longitudinal 

10 disturbance propagating along z, a periodic stress pulse 
propagating along y, and a dc response. 

The variables A^{k,s), B^(k,s), D^(k,s) etc. are 
potential constants which are determined by imposing the 
necessary boundary conditions and physical realizability 

15 requirements. For the system of Fig. 7 with H infinite, 
the terms involving C^(k,s) and (k,s} represent 
physically unreasonable solutions. For this reason, they 
are set to zero. The other six constants are determined 
with the boundary conditions. As stated earlier, these 

20 are that u and the normal components of the stress are 
continuous at y=0 and that the normal components of the 
stress vanish for y=-h. The stress tensor is calculated 
using the Dvihamel -Neumann relation shown in eq. (15) in 
order to properly take into account the contribution from 

25 the temperature grating. 



" """^[^ ] ^ f^""^ I - 2c44 )v . u - yrjSij 



(15) 



Using equations (8), (15) and (12)-(14), the 
boundary conditions can be imposed. The six relations 
thus derived are most conveniently expressed in matrix 
30 form, 
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£D=F 



(16) 



C = 



2lne-"^ -aine""* 

2in -2in 

-1 -1 

in -in 



lipe-P" 
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(21) 



m«, 1 



(22) 



r = Jl- 



(23) 



Pf^Tf 



(24) 



5 and and are the densities in the film and substrate 
respectively. The first two rows of the matrix equation 
ensure that the yy and the yz components of the stress 
respectively are zero at the surface of the film. The 
next two rows require continuity of the yy and yz 

10 components of the stress respectively at the interface. 
Finally, the last two rows require that the z and y 
components of the material velocities respectively are 
continuous at the interface. For notational convenience, 
the k and s dependence of the potential constants has 

15 been dropped. Also, the Laplace variable s has been 

written as s^--6)^=-v^K^ where w is the real frequency and 
V is the velocity of the acoustic mode solutions. This 
transformation limits the possible solutions to strictly 
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pseudO'Rayleigh modes with velocities below the 
transverse velocity of the substrate. A continuum of 
"leaky mode" solutions with complex frequency are also 
possible for velocities greater than v^g. These solutions 
are deunped as energy is "lost" to the semi-infinite 
substrate. Normal mode solutions of eg. (16) occur when 



Aside from kh and the matrix c only contains constants 
related to the film and the substrate. All of the laser 

10 and optical properties appear in the column matrix F. 
Thus, eg. (25) can be solved to yield the dispersion 
relations (gh vs. v where g is the acoustic wavevector) 
for the various pseudo-Rayleigh modes irrespective of the 
XSTS excitation source. 

15 An example of typical dispersion curves for the 

substrate loading case (ie. *Ts >^Tf) is shown for the 
lowest 16 modes in Fig. 8. As qh approaches zero only 
one mode, whose velocity approaches only one mode whose 
velocity approaches that of the substrate pure Rayleigh 

20 mode ^Rs, exists. At large gh, this lowest mode 

approaches the pure Rayleigh surface mode of the film 
with velocity '^Rf . Each of the higher velocity modes, 
often called Sezawa modes, has a limiting qh value below 
it does not propagate. At this qh value, the velocity of 

25 each mode is equal to the substrate transverse velocity 
^Ts. As shown in Fig. 8, as gh is increased, these modes 
all asymptotically approach the film transverse velocity. 

Xn order to solve for the displacements of a 
particular mode j excited through ISTS, one must solve 

30 the full eq- (16) for the potential constants at the qh 
value and velocity of interest. It can be shown using 
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eqs. (8) and (13-14) that the normal mode displacements 
with wavevector g ol>ey the following eguation 



uj(q.y.z.O = [ui(q.y)y + ui(q.y)i]c'<»'e'^'''* 



(26) 



qnj[Aie'"'"i -Bie-'«'"i]-. iq[Aj,c'"«'i -h < y < 0 



ui(q.y) = ' 



-qmjD;e-'*""i+iqDj,e-^i 



y>0 



(27) 
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ui(q.y) = 



-h<y<0 



(28) 



In these expressions, the tildes denote the residue of 
the quantity evaluated at the appropriate pole. Note 
that the material displacements oscillate with frequency 
5 Oj^v^q. Examples of the pseudo-Ray leigh displacements for 
various modes at one qh value are showing in Fig* 8. The 
total displacements at a particular wavevector and sample 
thiclcness will be the sum of the uj values for each 
pseudo-Rayleigh mode plus a nonoscillatory term due to 

10 steady-state thermal expansion. 
B- Probe Diffraction 

In this section, an analysis is presented of 
diffraction appropriate for the ISTS experiments 
described here in which ripple effects were found to ■ 

15 dominate. The equations below are the reflected and 
transmitted first order diffracted components of light 
from a corrugated half space, modified to deal with the 
case of arbitrary indices of refraction on either side of 
the corrugated interface, and to the case of standing 

20 wave corrugation. In addition, index modulation 



<WO 930 1 476A 1_l_> 



wo 93/01476 



PCr/US92/05679 



- 37 - 

contributions to the phase are neglected. The results 
are. 



.T- ic- .i«.k.R,r ,( {\'*'r(BO)coi^^ -(\'t(Bi)M; _ q(l r(90)(sin9, ^ sinS;) ) 
E, = -avIoC 3 . /,(a,)j^ " In^kj 0056^(0,58, * cos 9^), 



(29a) 



ai - S(i)niki(co$e^ +00585) 



(29b) 



20 



(29C) 



(l + r(ei))il2.co$e» +(l-r(ei))cosei q(l-»^r(0i)/ii2.5ine, +sin9il 
»i 5 \M I 

2nilCiCOs6j^-^cose| -cosQi j 



2cos8| 



(30a) 
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a', = 5(i)n,ki|^^cose^ - cosOi j ( 3 Ob) 



sine.+i = -Jl^sinei ± (30c) 



In these equations, 1^ is the intensity of the probe 
measured in medium 1, and is the wavevector of the 
5 probe measured in vacuum. 6^^ is the angle of the probe 
relative to the corrugated interface, a is the unit 
vector representing the polarization of the probe, r is 
the reflection coefficient for transmission from medium 1 
to medixim 2, and is the distance from the point at 

10 which the diffraction is generated to the detector* 

Also, g is the wavevector of the corrugation, and S{t) is 
the time dependent amplitude of the corrugation. 
Finally, the upper signs refer to the +1 diffracted 
order, while the lower signs refer to -1 diffracted 

15 order . 

For a thin supported film, there are four sources 
of first order ripple diffraction which must be 
considered. These are shown schematically in Pig. 9. 
Parts A and B of this figure represent the probe and 
20 diffracted beam paths for the two final signal directions 
seen experimentally. The first three paths show the 
components involving diffraction from the film surface 
ripple. The foxirth path describes the component arising 
from diffraction from the interface ripple. In general. 
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all four components contribute to the diffracted signal 
in the two signal directions. 



with these four diffraction contributions can be obtained 
5 using eqs. (29-30) , Snell's Law, and taking into account 
reflective losses at the film-air interface. In the 
following, it is .assumed that all light striking the 
film-substrate interface is reflected. In addition, due 
to the spatial size of the probe pulse, the thickness of 

10 the film, the index of refraction of the film, the 
wavevector of the corrugation, and the angle of 
incidence, it is a good approximation to assume that all 
of the diffracted beams give in Fig. 9a or 9b emerge from 
the same point on the film and travel parallel to one 

15 another. The results, written in terms of sums over 
pseudo-Ray leigh mode amplitudes, are shown in the 
following equations. 



The electric fields at the detector associated 



El = n^Kic'^e^'^i^^^ IuJy(q.y » -h)sin(vjqi) 



i-0 



(3la) 



(31b) 




(31c) 



i 



-ikiRog-inkihj £ujy(q,y -h)sin(v|qt) 
j=0 



(32a) 
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K sk fncosB coje J O » f^^»^)"'^o^^t <^ ' ^Ce.))cosQ, ^ q(l * f(9.))(nsin9, ^ sine,) ! (32b) 
^" *^ '\ 4cosei ' 4k,cose,{nco$0i -cosdi) J 



n 



(32c) 



h2- 



2h 



cos 9,* J 



(32d) 



E3 =+ai(ei)Vl7K3C^^e-"^i^oc-^i'»3 Zujy(q.y = -h)sw(Viqt) 

i=0 



(33a) 



(33b) 



sine'iiissnsinerfrt-^ 



(33c) 



h3 = 



2h 



cosS^ 



(33d) 
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erf^ = sin^^^isinej j (33e) 



E4 =Tat(ei)r(e*d)Vl^K4C^«c-*i'^oe-i'^«*'4 Iuiy(q.y = 0)sm(Vjqi) | . 



(34b) 



sine'd±i=sinerfri';^ 



(34c) 



Here, j runs from 0 to N and represents the mode 
solutions for j>0 and the dc term for j«0 (ie. Vo==0) . 
Also, t is the transmission coefficient for passage from 
the air into the film, and t' is the transmission 
10 coefficient for passage from the film into the air. 
Finally, in these expressions the small -argument 
expansion of the Bessel function (Jiix) x/2 for x«l) 
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has been used since the rippled amplitude (as estimated 
in section IV) times the grating vavevector magnitude is 
small compared to one for all excitation geometries 
investigated in this experiment. 
5 The signal intensity can be written as the square 

of the sum of the four electric fields given in egs. (31- 
34) . The resulting intensity will exhibit varying 
interference effects depending on the probing angle and 
wavelength, the index of refraction and thickness of the 

10 film, and the grating excitation wavevector. For the 
purpose of analyzing the relative intensities of the 
various pseudo-Ray leigh modes in an ISTS experiment, the 
only important interference occurs between the E-f ields 
derived from the surface ripple (£^-£3) and the field 

15 derived from the interface ripple (E4) . Such 

interference can cause a mode with a large surface 
corrugation amplitude to only appear as a weak Fourier 
component in the experimentally measured diffraction 
signal and vice-versa. The final ISTS intensity (I) can 

20 be written in a manner which emphasizes this effect as 

I = Io^l[{F(n.h.ki.q.ei)ujyCq,y = -h) + G(n.h.ki.q,ei)uJy(q.y = 0))sin(vjqo]x cc.J pgj 

where F and G are functions determined by summing terms 
in equations (31-34) . This expression was obtained by 
averaging over an optical cycle with the assumption that 

25 acoustic terms are constant over this interval. It can 
be seen that, in general, the measured signal consists of 
a dc term, a series of heterodyned components oscillating 
at the pseudo-Rayleigh mode frequencies which arise from 
a beating against the dc term, and terms oscillating at 

30 sums and differences of the normal mode frequencies. 
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C; Qualitative Considerations 

The relative intensities of the various modes seen 
with ISTS depend on two factors - the extent to which the 
different modes are excited by the pump pulses and the 
5 efficiency with which the surface and interface mode 
displacements diffract the probe pulse. Both of these 
factors can be quantified at a particular qh value by 
solving the full matrix eq. (16) for each mode. This was 
done using the elastic constants found for the 

10 Pyralin/ silicon system as described below. Aside from an 
arbitrary laser- intensity dependent constant factor, the 
only additional input parameter needed to perform this 
computation was the Pyralin absorption coefficient (• 
This was measured to be 1.3/Ltm-^ using a Pyralin film 

15 coated onto a fused silica substrate* Fig. 10 shows the 
results of the calculations at qh=2.5. This figure shows 
material displacements for the first eight modes that 
exist at this value of qh. The relative excitation 
efficiencies for the different modes are indicated by the 

20 relative amplitudes A. (A = 1 for the mode with the 
largest displacements.) The surface ripple amplitude 
relative to the maximum amplitude for each mode is given 
by R. Modes with large values of the produce AR 
contribute most to XSTS signal. There are several 

25 guidelines that can be determined through a study of 
these and other figures that allow for an understanding 
of the level of excitation and the resulting diffraction 
efficiency of a particular mode. First, modes with fewer 
spatial modulations along the y direction are excited 

30 more strongly than those with many nodes along y. This 
is reasonable since the heating excitation mechanism will 
tend to displace material in only one direction for a 
give value of z. Next, lower velocity modes in general 
have larger peak displacements than high velocity modes. 

35 This is due to the fact that as the mode velocity is 
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increased and approaches the substrate transverse 
velocity, the mode displacements become spread out over 
both the film and far into the substrate. This leads to 
small pesOc displacement amplitudes compared to those of 
5 the low-velocity modes whose displacements are well 

localized in the film. Finally, it is evident that some 
modes involve much more surface and interface ripple than 
others. Since detection results predominantly from 
diffraction from surface and interface ripple, only modes 

10 with significant surface displacements along y are 
detected. With these guidelines, it is possible to 
qualitatively explain relative signal amplitudes from the 
various modes over a wide range of qh. For purposes of 
data analysis, the calculations provide a clear 

15 indication of which modes in a thin film are likely to be 
observed in ISTS data. In most of the data, only one or 
two frequency components are observed predominantly* 
Dispersion curves like those in Fig* 8 can be fit to the 
frequency values determined from the data, which the 

20 elastic moduli as adjustable parameters to be deduced 
from the best fit. To do this correctly it is essential 
to know which modes are under observation, i.e. which of 
the many dispersion curves should be fit to the 
experimentally determined frequency values. Another 

25 qualitative consideration that explains why ripple 

amplitude is small for modes with large velocities (which 
explains drop off for modes at a given qh and for all 
modes as qh is decreased to small values) is that modes 
are fast because they involve substantial motion in the 

30 substrate which is very stiff compared to the film (i.e. 
it is motion in the substrate that "speeds the mode 10") . 
Since the substrate is very stiff, equal energy 
deposition will result in smaller displacements for fast 
modes . 
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Fig. 10 illustrates bulk and surface ripple 
amplitudes for a single value of wavevector q, given the 
film thickness h. Fig. 11 is a three-dimensional plot 
showing the calculated film surface ripple (the product 
5 AR) excited through ISTS as a function of qh and velocity 
for the various modes in 1 fivi film. Fig. 12 displays the 
analogous results for the interface ripple. The x-y 
plane of these figures shows the mode velocity dispersion 
as in Fig. 8, while the z axis shows the relative ripple 

10 displacements for a given energy density of the ISTS 
excitation pulses. Note that since the displacements 
associated with each pseudo-Ray leigh mode are different, 
the relative ripple displacements for the film surface 
and the film * substrate interface are also different. 

15 Both, however, exhibit the same general trend of 

decreasing ripple amplitude with increasing qh. For 
fixed h as in this simulation, this corresponds to 
decreasing mode displacement amplitude with increasing 
wavevector which, as discussed earlier, is expected for 

20 ISTS excitation. In addition to this decrease in 

amplitude which is approximately like l/q for values of 
qh above 1, there is a distinct decrease in the magnitude 
of the surfaces ripple for small values of qh. This 
decrease is attributed to the fact that as q is decreased 

25 to values near zero, even the lowest modes take on large 
velocities. Thus, as discussed above, these modes are 
not well localized and so have small peak displacements. 
For this reason, the ripple due to such modes is small. ^ 
In addition to the decrease in surface ripple at 

30 both high and very low values of q, it is seen in Fig. 11 
that there is a "cross-over" near qh=1.5 where, upon 
increase in qh, the surface ripple due to mode two 
becomes small and that due to mode one becomes large. 
This occurs because the character of these two m des 

3 5 switch in this region. This "crossover" is evident from 
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the dispersion curves of Fig. 14, in which there is an 
"avoided crossing" of the first two frequencies, and 
through a comparison of modes one and two at gh=2 • 5 given 
in Fig. 10 and at gh^O.S given in Fig. 13- (The presence 
5 of this cross over effect is a property of ripple 
phenomena. } 

The relative contribution of each pseudo-Ray leigh 
mode to ISTS signal can be calculated using eq. (16) and 
the ripple information contained in Figs. (11-12). The 

10 resulting mode intensity ratios can then be compared to 
the ratios found in the raw ISTS data for each value of 
qh and h. When the interface displacements are not 
negligibly small compared to the surface displacements, 
interference effects between these two contributions to 

15 diffracted signal become important so that with different 
experimental conditions - such as probe angle of 
incidence or film thickness - ISTS data can show vastly 
different mode intensity ratios (cf . eq, 35) . 
Conversely, when one or the other set of displacements 

20 dominates, interference effects can be ignored and the 

resulting mode intensity ratios become independent of the 
precise experimental conditions and mimic the behavior of 
the relevant ripple displacements. From Figs. (11-12) it 
is clear that the theoretical surface displacements are 

25 over an order of magnitude larger than the interface 

displacements for the two lowest velocity modes with qh 
values above -0.5, while for the other modes and at lower 
qh, the surface and interface displacements are of the 
same order. Calculations performed at other film 

30 thicknesses between 1 and 5 show this same general 
behavior. Thus, for the samples examined in this study, 
one would expect the relative intensities of the two 
lowest modes to mirror the relative intensities shown in 
figure 9 for qh>0.5. (An analysis of the complete mode 
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intensity spectrum requires the full solution of eq. 
(35).) 

D. Quantitative Analysis of ISTS Data 

As stated earlier, the dc term dominates the ISTS 
5 signal for these experiments. In fact the power in the 
y zero-frequency peak of the Fourier transforms ranged from 

30 to 200 times greater than the acoustic peaks. Thus, 
the heterodyne terms in equation (35) dominate and the 
frec[uencies of the most intense Fourier peaks will 

10 correspond to the fundamental pseudo-I^ayleigh mode 
frequencies Oj. 

The values of found for all the scattering 
wavevectors and sample thicknesses in this study were 
converted to phase velocities and plotted versus qh. The 

15 results are showing in Fig. 14. in cases where multiple 
modes were excited^ only those with Fourier transformed 
intensities greater than half that of the largest were 
plotted. In addition, to eliminate the possibility of 
spurious peaks due to sample imperfections, data points 

20 were taken at several spots on each sample and only 

frequencies that were consistent from spot to spot on the 
same were kept. The estimated experimental uncertainties 
for the velocities range from 1% at the highest qh values 
to 5% at the lowest qh values observed for each sample. 

25 The qh values are accurate to ±0.1. The main source of 
error results from uncertainties in the scattering 
wavevector which are derived from uncertainties in the 
mechanical measurement of the scattering angle. This 
error can be substantially reduced by deducing the 

30 scattering wavevector directly from ISTS data taken with 
a well characterized reference sample. 

The solution for the pseudo-Ray leigh mode 
* dispersion for this system requires the density and the 

longitudinal and transverse velocities for both the 

35 silicon and the Pyralin as inputs. For silicon, pg= 2.33 
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g/cm and the isotr pic velocities are ^Lg* 8945in/s, and 
^Ta=5341 m/s. For Pyralin, the values given by Dupont for 
the density and Young's modulus (Y) are 1.45 g/cm^ and 
2.4 GPa respectively.^^ To the inventor's knowledge, 
there is no reported measurement of Poisson's ration (v) 
for this film. All three of these values are necessary 
to obtain the longitudinal and transverse velocities for 
the film according to egs. (36) and (37). 



v2 = 



^~2p(l+v) " (36) 



10 „2 _ 

P(l + V) 



Different sets of dispersion curves were generated by 
solving eg. (25) using the silicon and Pyralin parameters 
given above and varying Poisson's ration for Pyralin from 
0 to 0.5. None of these theoretical data sets were able 

15 to fit the experimental data adequately. The 

experimental data was thus fit numerically using the 
Marquardt-Levenberg nonlinear least squares algorithm and 
allowing all three of the Pyralin parameters to vary. In 
performing such a calculation, one must assign each 

20 experimental point to a particular dispersion curve. 
This was fairly straightforward for the data shown in 
Fig. 14 since the experimental points group together into 
well-defined curves resembling the pseudo-Ray leigh 
dispersion curves. The theoretical results shown in 

25 Figs. 11 and 12 and discussed above indicate that the 

lowest velocity pseudo-Rauleigh modes should give rise to 
the strongest signals for most of the wavevectors used. 
It was therefore assumed that the two lowest-velocity 
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sets of data points correspond to the two lowest-velocity 
pseudo-Ray leigh dispersion curves. The Pyralin 
velocities (^Lf and ^Tf) parameters were adjusted to fit 
the points on these lowest-velocity curves with the 
5 density fixed at the measured value of 1.45 g/cm^. Fig. 

* 14 shows the data points and the fit. Using the values 
obtained, the other pseudo-Ray leigh dispersion curves 
were calculated and plotted on Fig- 14. The agreement 
between these higher-velocity pseudo-Ray leigh dispersion 

10 curves and the higher-velocity data points is excellent, 
with no fvirther adjusting of any parameters* 

The Pyralin velocity parameters determined from 
the fit are ""Lf » 2650 m/sec and ""Tf^ 1160 m/sec. While 
it is never certain that the global minimum in the value 

15 of using nonlinear least squares fitting routines is 
reached, these values correspond to the lowest minimum 
found after starting from many different initial guesses 
for the parameter values. Varying either of the Pyralin 
velocity parameters by 10% led to substantially worse 

20 fits. Varying the substrate parameters by the same 

amount led to worse fits for values of gh<l, but did not 
affect the quality of the fits at higher qh values. It 
is thus estimated that the overall uncertainties in the 
pyrlain velocity parameter values to be ±5%. 

25 The relative intensities of the two lowest modes 

can be compared to the film surface displacements of Fig. 
11 for qh>0.5. Experimentally, Fig. 14 shows that at 
least one of these two modes is present for all the 
samples. According to figure 9, the second lowest mode 

30 should dominate from qh — 0.5 to qh-1.5 while at higher qh 
values, i.e. above the crossover discussed earlier, the 
lowest mode should have a higher intensity. Upon 

* examination of figure 12, this qualitative behavior 
definitely obtains since, for all the samples, the 

35 density of observed data points is much higher for the 
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second mode 'than for 'the first gh values less than the 
crossover, while at higher gh values the only observed 
dalia points correspond to the first node. A cjuantitative 
mode intensity analysis is possible by accura-tely fixing 
5 or measure the probe smgle of incidence 6^ since, 

according to eqs (29-34), changes in 6^ significantly 
change the values of F and G in eg. (35) which in turn 
can change the final mode intensities seen in an ISTS 
experiment. 

10 The outlying high velocity data points observed 

are not strictly pseudo-Rayleigh surface modes but rather 
represent modes which arise from interactions with the 
free substirate boundary. Their counterparts for the 
semi-inf init:e substratie are known as "leaky" modes which 

15 are damped due to ^e fact that they lose energy to the 
semi-inf init.e substrate as they propagate. The behavior 
of the modes can be elucidated by treating the finite 
size of the substrate explicitly. This is demonstrated 
below. 

20 The solution for the transient grating excitation 

of a thin film coating is generalized to include the 
effects of a substrate with finite thickness. The formal 
solution for the displacement potentials is still given 
by eqs. (13-14). However, now the finite substrate has a 

25 free boundary at y=H (cf . Fig. 7) which is ignored for 
the semi-infinite system. For this reason, the and 
terms of eg. (13-14) eure no longer unphysical and must be 
retained in order that the normal components of the 
stress at this new surface can be fixed to zero. 

30 With the two new boundary condition equations and 

the two new unknown potential constants, eq. (16) becomes 
an 8x8 matrix equation. Using the Duhamel-Neiunann 
expression for the stress (eq.) 15), and assuming as 
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before that there is no heating of the silicon, it can be 
shown that the new matrix equation takes the form 

fiV=P* (38) 

where 
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-(6(k).6(lc-q)) vt,^f -(vu^-2vTr^) .(;f ^ , ^ 



8(k) ^ scc-q) 



? 



f-(vui-2vTT^ cM S(k) , 8(k-q) 



Yxr^pkv 





5(k) 


ik^ 






S(k) 








5(k) 







6(k-q) 



8(lc-<l) 



a(k-q) 



(41) 
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Information is extracted from this matrix equation 
with the same techniques that were used for the 6x6 
equation in the main text. The resulting dispersion 
curves for the thirty lowest velocity modes of a system 
5 composed of 3 fim Pyralin film on a 330/im silicon 
substrate are showing in Fig. Al. One can see that 
relaxing the semi-infinite substrate assumption give rise 
to extra plate modes propagating predominantly within the 
substrate with velocities above the substrate transverse 
10 velocity ^T^ which was, in the semi-infinite case, the 
cutoff for propagating modes of the system. The 
dispersion curves for velocities sufficiently below the 
'^Tg are almost identical to those for the pseudo-Ray leigh 
modes calculated using the semi -infinite substrate system 
15 and shown in Fig. 14. The only differences occur at 
phase velocities near the cut-off and, for the lowest 
velocity mode, at very low (<0.2) qh values. This 
behavior is expected and confirms that the semi-infinite 
substrate approximation is adequate to explain nearly all 
20 of the pseudo-Ray leigh modes seen in these experiments. 

The extra plate modes above ^Tg are produced as a 
result of acoustic reflections from the free substrate 
s\irface. They can be obtained within the semi-infinite 
substrate approximation as well, but only as "leaky" 
25 modes with finite lifetime due to energy flow into the 

infinite substrate. The existence of these modes for the 
Pyralin/ silicon system used in these experiments provides 
an explanation for the ISTS data points occurring above 
'^Tf in Fig. 14. 
30 £. Discussion 

The values of ^Lf = 2650+130 m/sec and ^Tf = 1160+ 
60 m/sec deduced in this study for Pyralin 2555 combined 
with the density of 1.45 g/cxa? correspond to a Young's 
modulus of 5.4+0.5 GPa and a Poisson's ratio of 
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0.38+0*02 . (The uncertainties given for Young's modulus 
and the Poisson ratio are calculated from the 
uncertainties for ^Lf and ^Tf determined from the study 
above and do not account for uncertainties in the film 
5 density.) To the inventor's knowledge, this represents 
the first measurement of Poisson 's ratio in this 
material. This value for Young's modulus is over 100% 
higher than the previously measured value of 2.4 GPa. 
The value of Young's modulus can depend very sensitively 

10 on the method of film curing. Thus, the difference is 
believed due to intrinsic difference in sample 
characteristics . 

Polyimide films can exhibit variation in the 
degree of chain orientation and in density as a function 

15 of depth. These measurements yield values which are 
averaged over the film thickness and so do not provide 
direct information about depth dependences. Measurements 
of thinner films may yield different parameters, 
indicative of differences in film properties very close 

20 to the substrate surface. By tilting the grating 

wavevector so that it has a component perpendicular to 
the film, it is possible to determine separate in-plane 
and out-of-plane elastic moduli. The theoretical 
treatment outlined here could be generalized to take this 

25 into account. 

In general, there are several different protocols 
available for extracting the elastic parameters from ISTS 
data on thin supported films. One method, which was 
demonstrated above, involves determining the mode 

30 velocities at a range of gh values - by either changing 
scattering angle or sample thickness or both - and then 
varying the elastic parameters until a good match between 
theoretical dispersion curves and data is obtained. 
Experimental error can be reduced with improved accuracy 

35 of the scattering wavevector measurement and an analysis 
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along these lines with fewer gh values and higher 
accuracy can be achieved. An alternative or 
complementaxy method involves using the relative ISTS 
signal intensity information obtained for the different 
5 modes at each gh values. As discussed, a guantitative 
analysis of the relative intensities reguires accurate 
specification of the incident probe angle, the indices of 
refraction of the film and substrate, and the film 
thickness. The first parameter is easily determined in 

10 the ISTS experiment. The indices of refraction and film 
thickness can be measured by other technigues, or treated 
as free parameters in a fitting scheme. The elastic 
parameters may be extracted using data from only one gh 
value at which multiple modes are observed, by fitting 

15 both the pseudo-Rayleigh mode velocities and relative 
intensities. With efficient computing algorithms to 
accomplish fitting, the elastic p8urameters may be 
extracted with the same real-time rates demonstrated for 
the data acguisition. Generally, the experiment can be 

20 optimized with an improved probe pulse temporal profile. 
£• Determination of Adhesion Oualitv 

In the analysis described above, the film and 
substrate were modeled such that together they form an 
acoustic waveguide and together determine the freguencies 

25 of the acoustic waveguide modes that propagate in the 
film-substrate system. Referring to Figs. 16 and 16a, 
comparing the dispersion curves for a tightly bound film- 
substrate system with a system in which the substrate is 
removed (which would be eguivalent to a film-substrate 

30 system with a substrate that has a very low stiffness) it 
is seen that the waveguide freguencies change 
dramatically. These changes can be used to detect 
regions where the film has become debound from the 
substrate by sampling the film at various locations. 
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Further, the degree of film substrate adhesion can 
be determined. For example, the absence of a substrate 
is equivalent to the case where ^e film is no longer 
stuck to the substrate. If the film is slightly stuck 
5 but not tightly bound, the acoustic waveguide frequencies 
are intermediate between the case where the substrate is 
present and where it is removed. By measuring the actual 
positions of the frequencies, the degree of adhesion can 
be determined. In the adhesion model, two parameters are 

10 introduced (k^* ky, further discussed below) . One is a 
spring stiffness parameter per unit area for motions 
parallel to the plane of the film (k^) . The other is a 
spring stiffness parameter per unit area for motions 
perpendicular to the film (ky) • Although each pareuneter 

15 affects the waveguide frequencies, for films with 

intermediate degrees of adhesion, we expect that since 
the film is in contact with the substrate, the stiffness 
parameter (ky) for motions perpendicular to the film will 
be large (i.e. close to the tightly bound case). 

20 Referring to Fig. 17 only the stiffness parameter for 
motions parallel to the plane (k,) of the film has been 
varied. (The plot wherein k^^ infinity corresponds to 
the tightly bound case. Fig. 16, and the case wherein 
k^^^O corresponds to the case of a film on a frictionless 

25 surface.) As can be seen from this figure, there is a 
smooth variation in the acoustic waveguide mode 
frequencies between the tightly bound case and the i 
frictionless interface case (i.e. the case where there is 
no resistance to film motions parallel to the plane of 

30 the film but there is resistance to motions perpendicular 
to the film) . The waveguide frequencies are most 
sensitive to changes in the adhesion at small wavevector 
times thickness values. Therefore, this region can be 
probed (by varying the wavevector of the excitation) . 

SUBSTITUTE SHEET 



BNSDOCID: <WO ^9301476A1J_> 



wo 93/01476 



PCT/US92/05679 



- 58 - 

The plo'ts in Figs. 16 and 17 were genera'ted by compu'ter 
algorithm using inputs as follows. 



(42) 



5^ (d^) * Cd") (42a) 

(42b) 
(42C) 

where is the stress tensor; Uj^ is the displacement 
vector; la^ is an elastic stiffness/unit area for motion 
along the z axis; ICy is an elastic stiff ness/unit area 

10 for motion along the y axis (perpendicular to the plane 
of the sample s\irface) . These boundary conditions are 
used with the equations of motion to determine acoustic 
mode dispersion and excitation efficiency. (Discussion 
of models for a different purpose can be found in F.J. 

15 Margetan et al.. Journal of Nondestructive Evaluation . 
Vol. 7, Nos. 3 and 4, 1988, p. 131.) 

An inertial mass term was used to describe excess 
mass at the crack is dropped for the adhesion model. 
F. Analysis of Thermal Diffusion Data 

20 Referring to Fig. 18, the following is a 

derivation of a technique for measurement of thermal 
diffusion, using boundary conditions for a thin film. 
Starting with the equation of motion: 
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cpdT 

K 3t 



Q 

K 



(43) 



Q represents the heat source; K is the thermal 
conductivity; C Is the heat capacity per unit mass; and p 
is the density. 
5 It is assumed that heat flow occurs within the 

film and out of the film into the air and substrate. The 
heat flow out of the film obeys Newton's laws of heat 
propagation (i,e, heat flow is proportional to the 
temperature gradient leaving the following boundary 
10 conditions on heat flow: 



ay 



y-0 



ar 
dy 



(44) 



In the above, the z axis is in the plane of film 
and the y axis parallel to the film thickness with the 
origin, y=0, of the film surface and where the film 
15 thickness is y=d| h and h' relate to the degree of 

thermal coupling between the film and air and film and 
substrate, respectively. For ISTS excitation: 

With Q=A(l+e^«*) (e~fy)y(t) , the transformed 
solution (using transforms above) is; 



T(s.k.y) = B(s.k)exp|-y[k2 ..^J| + C(s.k)exp|y(^k^ +^J| 

SUBSTITUTE SHEET 



(45) 



BNSDOCID: <WO 9301478A1_L> 



wo 93/01476 




PCr/US92/05679 



- 60 - 



S(k)-fS(k-Jp) 

K J 



(46) 



where A Is tiie ampli'kude factor related to excitation 
pulse intensity and optical absorption coefficient, q is 
the wave vector and f is the optical absorption 
coefficient. Applying the boundary conditions and 
inverting the transform yields 



T(y.2.t)=j^l 







r 

lexp 


[ CP J 







cpd' 



r„) 



(47) 



with 



(48) 



and 



2A 










2-d(h*-h)smr6 +(h'-h)j^^-dcotrB jcosfaj 



(49) 
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In the quasi-static limit, the displacements are 
related to the temperature as follows: 




(50) 



5 where is the coefficient of thermal expansion and Jl is 
the displacement. 

In this limit, the temporal dependence of Jl is 
determined by the temporal dependence of T. Thus, the 
signal is a sum of exponentials with the same decay rate 
10 in equation (46) • The different decay rates can be 
extracted with a linear prediction routine. With the 
rates determined in this manner, comparison can be made 
to theory to determine: 

C( (51) 



15 by changing q and d, then, predictions can be made for 
the form of the experimental decay. In this way, 
equation (51) can be compared to experimental values to 
determine the accuracy of the calculated values. 
Further Embodiments 

20 It will be evident from the above that the 

invention enables many embodiments and advantages. While 
the experiments may be carried out on thin films to 
particular advantage as discussed above where surface 
ripple is the predominant phenomenon leading to 
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diffraction of the probe, the ripple effect may also be 
induced on thicker samples. The experiment can be 
optimized to selectively analyze diffraction from ripple 
over diffraction from bulk, refractive index 
5 modification. The ripple effect can be preferentially 
detected by maximizing the amount of reflected light by 
increasing the incident angle of the probe beam. The 
optimum incident angle may be determined based on sample 
thickness and the refractive index to enable a large 

10 incident angle while avoiding substantial losses due to 
total internal refraction. The detection of diffraction 
from the ripple effect can also be optimized on 
freestanding samples by detection of transmitted and 
reflected probe radiation and subtracting the former, 

15 which contains bulk-properties information, from the 

latter, which contains ripple information. (It will also 
be understood that for thin films, ripple from both sides 
of a film can be detected, even for films on a substrate, 
in which case ripple of the sample-substrate interface 

20 can be detected.) The presence of signal induced from 
bulk properties can also be detected by study of the 
polarization of the diffracted beam. A polarizer can be 
positioned after the probe focusing optic (e.g. optic 54) 
and another polarizer positioned before the detector 

25 (e.g. detector 58). As discussed above bulk properties 
vary the incident polarization. The polarization can be 
studied as a function of incident angle to determine the 
optimum angle for detecting ripple induced diffraction. 
Further, the excitation radiation can be optimized to 

30 enhance the diffraction signal arising from ripple. The 
wavevector dependence of the signal can be studied (e.g. 
with a computer) to determine that the diffraction signal 
is predominately from ripple. Generally, a smaller 
excitation wavevector (longer excitation wavelength or 

35 smaller angle of incidence of the excitation beams) 
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enhances the ripple effect. In the experiments on the 
sample as described above, the maximum vavevector for 
which signal was obtained was about q=1.8lAtm-^. Larger 
maximum wavevectors can be used if a less stiff sample or 
higher intensity probe laser is used. 

other embodiments are within the scope of the 
following claims. 

What is claimed is: 



SUBSTITUTE SHEET 



BNSDOCID: <WO_9301476A1J_> 



wo 93/01476 




PCT/L'S92/05679 



- 64 - 
Claims 

1. An apparatus for measuring the properties of a 
sample of material, comprising: 

a first, excitation source for producing 
5 excitation radiation adapted to impinge upon said sample 
of material, said excitation radiation comprising pulsed 
radiation composed of at least two component pulses which 
interfere within said sample, said excitation radiation 
being sufficient to induce a transient phonon in said 
10 material which gives rise to a transient, time dependent 
periodic ripple morphology on a surface of said sample, 

a second, probe source for producing probe 
radiation arranged to reflect from the periodic ripple 
morphology on the surface of said sample to form a 
15 diffraction signal, 

a detector for detecting the diffraction signal 
from said probe source radiation reflected from said 
surface , and 

an analyzer for selectively analyzing said 
20 diffraction signal formed by said transient ripple 
morphology. 



2. The apparatus of claim 1 wherein said probe 
source produces radiation that is not substantially 
absorbed by said sample. 

25 3 . The apparatus of claim 1 wherein said probe 

source produces radiation that is absorbed about 10% or 
less than said radiation from said excitation source. 



4. The apparatus of claim 3 wherein said probe 
source produces radiation that is absorbed about 1% or 
30 less than said radiation from said excitation source. 
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5* The apparatus of claim 1 wherein said 
excitation source radiation comprises ultraviolet 
radiation absorbed sufficiently by said sample to induce 
heating that gives rise to said phonon and ripple 
5 morphology • 

6. The apparatus of claim 5 wherein said probe 
source radiation comprises visible radiation not 
substantially absorbed by said sample. 

7. The apparatus of claim 1 wherein said sample 
10 comprises a thin sample of about 500 fm or less. 

8. The apparatus of claim 7 wherein said thin 
sample is about 10 fim or less. 

9. The apparatus of claim 1 wherein said sample 
comprises a thin polymeric film. 

15 10. The apparatus of claim 9 wherein said film 

comprises a free-rstanding film. 

11. The apparatus of claim 9 where said film is 
disposed on a support. 

12. The apparatus of claim 1 wherein said 

20 analyzer includes a polarizer for determining change of 
polarization of said probe beam after diffraction from 
said surface to selectively analyze diffraction from said 
ripple morphology. 

13. The apparatus of claim 1 wherein said 

25 analyzer is for analyzing said signal as a function of 
the wavevector and to selectively analyze diffraction 
signal from said ripple morphology. 
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14, The apparatus of claim l further Including a 
probe beam and detector constructed for detecting and 
resolving a substantial time-portion of the time- 
dependent diffraction induced by each excitation pulse. 

5 15. The apparatus of claim 14 wherein said probe 

radiation is of a selected pulse width and said detector 
is adapted to detect the time-dependent diffraction for 
the duration of said probe pulse. 

16. The apparatus of claim 15 wherein said 

10 detector is constructed to detect the entire detectable 
time-dependent diffraction induced by each excitation 
pulse • 

17. The appauratus of claim 16 wherein said 
excitation radiation has a pulse width on the order of 

15 psec duration and the probe radiation pulse width is on 
the order of nsec duration. 

18 • The apparatus of claim 17 wherein said 
detector has a time resolution on the order of 1 nsec. 

19. The apparatus of claim 14 wherein said probe 
20 pulse has a peak power output of about 1000 watts or 

greater . 

20. The apparatus of claim 19 wherein said laser 
comprises a Q-switch YAG laser. 

21. The apparatus of claim 14 wherein said 
25 excitation pulse is generated from a pulsed laser and 

said probe pulse is generated from a CW laser. 
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22. A method for measuring the properties of a 
sample comprising: 

impinging a pulse of excitation radiation upon 
said sample, said excitation radiation being composed of 
5 at least two component pulses which interfere within said 
sample, and being selected to induce a transient phonon 
in said sample which gives rise to a transient, time 
dependent ripple morphology on a surface of said sample, 

reflecting probe radiation from the periodic 
10 ripple morphology on the surface of said sample to form a 
diffraction signal, and 

selectively detecting the diffraction signal from 
said probe source radiation reflected from said ripple 
morphology* 

15 23. The method of claim 22 wherein said sample is 

• about 500 /im thick or less. 

24. The method of claim 23 wherein said sample is 
about 10 ^m thick or less. 

25. The method of claim 24 wherein said sample is 
20 a pure polymer sample. 

26. The method of claim 25 comprising selecting 
excitation radiation that is absorbed by the sample. 

27. The method of claim 26 comprising selecting 
probe radiation that is absorbed less than said 

25 excitation radiation. 

28. The method of claim 27 comprising selecting 
probe radiation that is not substantially absorbed by 
said sample. 

SUBSTITUTE SHEET 



QNStX)CID: <WO__S301476A1J_> 



wo 93/01476 




PCr/US92/05679 



- 68 - 

29. The method of claim 28 wherein said 
excitation radiation comprises UV radiation and said 
probe radiation comprises visible radiation. 



30. The method of claim 22 comprising detecting 
5 and resolving a substantial time-portion of the time* 

dependent diffraction induced by each excitation pulse. 

31. The method of claim 30 wherein said probe 
source is a CW laser. 



32. The method of claim 30 further comprising: 
10 providing a probe source having a power output of 

about 1000 watts or greater, and 

detecting and resolving a substantial time-portion 
of the time-dependent diffraction induced by each 
excitation pulse. 

15 33. The method of claim 30 further including 

signal averaging the diffracted radiation from multiple 
excitation pluses. 

34. The method of claim 22 further comprising 
analyzing said signal to selectively analyze diffraction 

20 from said ripple. 

35. The method of claim 32 comprising analyzing 
the polarization of said diffracted radiation. 

36. The method of claim 34 comprising analyzing 
the diffracted radiation as a function of wavevector. 

25 37. The method of any one of claims 33 to 35 

including optimizing the diffraction from said ripple by 
varying the angle of incidence of said probe beam. 
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38. The method of any one of claims 33 to 35 
including optimizing the diffraction from said ripple by 
varying the wavelength or angle of incidence of said 
excitation radiation. 

5 39, An apparatus for measuring the properties of 

a thin sample of polymeric material comprising: 

a first, excitation laser source for producing a 
pulse of radiation adapted to impinge upon said sample, 
said excitation source comprising a pulsed source 
10 composed of at least two component pulses which interfere 
within said thin film, said excitation radiation adapted 
to induce a transient phonon in said material, 

a second, probe laser source for producing 
radiation, said probe radiation being of selected 
15 wavelength not substantially absorbed by said sample and 
arranged to reflect from the surface of said sample to 
form a diffraction signal, and 

a detector for detecting the diffraction signal. 

40. The apparatus of claim 39 further including a 
20 probe beam and detector constructed for detecting and 

resolving a substantial time*portion of the time- 
dependent diffraction induced by each excitation pulse. 

41. The apparatus of claim 40 further including a 
signal average for signal averaging signal from multiple 

25 excitation pulses. 

42. The apparatus of claim 40 wherein said 
excitation pulse is on the order of psec duration and the 
probe pulse on the order of nsec duration. 

43. The apparatus of claim 42 wherein the probe 
30 laser source peak power is about 1000 watts or more. 
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44* The apparatus of claim 43 wherein the peak 
power is about 10,000 watts. 

45. The apparatus of claim 43 wherein said laser 
is a YAG Q-switched laser. 

5 46. The apparatus of claim 45 wherein said 

detector has a response time on the order of 1 nsec. 

47. The apparatus of claim 39 wherein said 
diffracted radiation is reflected from the back surface 
of said sample opposite said radiation sources. 

10 48. The apparatus of claim 39 wherein said probe 

source is a CW laser. 

49. The apparatus of claim 48 wherein said probe 
source is a gated argon ion laser. 

50. The apparatus of any one of claims 39 or 48 
15 wherein said probe has a power output of around 1 watt. 

51. The apparatus of claim 39 including an 
analyzer adapted for determining the adhesion of said 
polymer sample on a substrate surface from said 
diffraction signal. 

20 52. The apparatus of claim 39 including an 

analyzer adapted for determining the thermal diffusion of 
said sample from said diffraction signal. 
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